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TEMPERATURE DEPENDENCE OF THE RATE CONSTANTS 
FOR THE REACTIONS 0 + + 0 2 - 0 2 + + O 
AND 0 + + N 2 - NO + + N 
Peter Stubbe* 

ABSTRACT 

The dependency of the rate constant k r for the reactions 

0 + + o 2 - o 2 + + 0 k x (1) 

0 + + N 2 - N0 + + N k 2 ( 2) 

on the gas temperature T, the vibrational temperature T v of the molecular 
reactant, and the relative flow velocity v s between the gases, is studied theo- 
retically by means of a collision theory approach. The resulting expression for 
k r contains the two unknown parameters A (activation energy) and P 0 (scaling 
factor) which have to be determined empirically. k x (T) is found to have a min- 
imum at about 650°Kand a nearly linear increase for temperatures above 750°K. 
For reaction (2), the rate constant is presented as a function of T and A, since 
uncertainties in the experimental results do not allow a reliable determination 
of A. The increase rate ofk 2 with T v (N 2 ) , the vibrational temperature ofN 2 , 
falls off strongly with increasing gas temperature. 

*NAS-NRC postdoctoral resident research associate on leave of absence from the Max-Planck- 
Institut fuer Aeronomie, Lindau/Harz, Germany. 
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TEMPERATURE DEPENDENCE OF THE RATE CONSTANTS 


FOR THE REACTIONS 0 + + 0 2 - 0 2 + + O 
AND 0 + + N 2 - NO + + N 

1. INTRODUCTION 

The chemical reactions of 0 + with 0 2 and N 2 , 

0 + + 0 2 - o 2 + + 0 k x (1) 

0 + + N 2 - N0 + + N k 2 , (2) 

provide the main electron-ion loss mechanism in the ionospheric F-region. A 
great number of laboratory measurements of the reaction rates k x and k 2 at a 
temperature of 300°K (summarized by Ferguson, 1967) yielded ionospheric loss 
coefficients which led to serious difficulties in explaining some observed F-layer 
phenomena, mainly the maintenance of the nighttime F-layer. Therefore, it 
would be highly desirable to find a negative temperature dependence of k x and k 2 
as claimed by Donahue (1966) and Hanson and Cohen (1968) . 

Recent measurements of k l and k 2 as a function of temperature between 
77°K and 600°K (Nakshbandi and Hasted, 1967; Smith and Fouracre, 1968; Dunkin, 
Fehsenfeld, Sehmeltekopf, and Ferguson, 1968) showed very consistently that k t 
is indeed decreasing with increasing temperature. For k 2 , on the other hand, 
the situation is not as clear. Nakashbandi and Hasted (1967) report a slight in- 
crease of k 2 with temperature, while after Dunkin et al. (1968) k 2 clearly shows 
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a negative temperature dependence. According to Warneck (1967), finally, k 2 
is constant between 700°K and 1400°K. 

The temperature dependence of and k 2 described so far applies to a sys- 
tem in thermal equilibrium, i.e. with a distribution of excited states as given by 
Boltzmann's law. According to Thomas and Norton (1966), k 2 may be appreciably 
higher when the vibrational states of N 2 are excessively excited. This effect was 
measured quantitatively by Schmeltekopf, Fehsenfeld, Gilman, and Ferguson 
(1967) with the result that for a gas temperature of 300°K k 2 strongly increases 
with the vibrational temperature T v . As pointed out by Thomas and Norton (1966), 
this process may be of importance for the ionosphere during magnetic storms. 

The purpose of this paper is to present a theoretical model which is capable 
of reproducing the above cited experimental results as well as making predictions 
for physical conditions and temperature ranges which are not covered by the ex- 
periments. In the light of recent ideas about large s<^ale horizontal neutral gas 
motions (King and Kohl, 1965) it appears that a considerable amount of kinetic 
energy is stored in the atmosphere. It seems worthwhile to study how this could 
influence the rate of an ion-molecule reaction or, in other words, how the reac- 
tion rate depends on the relative velocity between the ion and neutral gas. We 
will include this question in our considerations. 

Our theory is based on a collision theory approach which is essentially a 
manifold modification of the Gioumousis-Stevenson theory (Gioumousis and 
Stevenson, 1958) that was successfully applied to some ion-molecule reactions 
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involving hydrogen. Additional information concerning the temperature depend- 
ence of the steric factor will be taken from the transition- state theory. 

2. THEORY 

A collision between two particles can lead to a chemical reaction when two 
conditions are fulfilled. First, the relative kinetic energy at the moment of the 
imp act must exceed the activation energy A. Second, the particles must be so 
oriented relative to each other that the reacting groups come close to each other 
during the collision. The fraction of collisions with a proper orientation of the 
particles is represented by the steric factor P. The determination of the chemical 
reaction rate k r by means of collision theory simply amounts to a count of the 
number of collisions that satisfy the above conditions. A theoretical prediction 
of the constant activation energy A and the temperature dependent steric factor 
P(T) lies beyond the scope of this theory. There is, however, another theoretical 
approach to the problem, provided by the transition-state theory. This theory 
enables a description of the temperature dependence of P. Thus, adding this 
information to our collision theory derivation of k r , we are left with only two 
unknown parameters, namely the activation energy A and the constant P Q in the 
steric factor P. A and P 0 must be determined empirically. 

According to these preliminary remarks, k r is defined by 

dN i dN n 

dk = PQ(g) g n N (3) 

in 
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where g is the relative velocity between an ion and a neutral particle. Q(g) is 
the effective cross section for collisions with relative kinetic energies greater 
than A. N. and N n are the numbers of ions and neutral particles per cm 3 , re- 
spectively. dN. dN n ^N. • N n is the number of ion-neutral pairs with a relative 
velocity between g and g + dg, related to the total number of ion-neutral pairs. 
3 n 7 dN n /N i • N n can be obtained by integration over 6 in Equation (47) of Stubbe 
(1968) as 


dN. dN 


UN. dIN /P\3/2 _ 2 _ ■ ■ i i 

N7n7 = ^in) e V * e_Rg g 2 d e (2n + l)T ( 2Rv s g) 2n 

n= 0 


(4) 


This expression is based on the assumption that both gases possess Maxwellian 
velocity distributions which are displaced by the macroscopic drift velocities v. 
and y n , respectively (see Stubbe, 1968, Equations (44a, b)). v g andR are defined 
by 


V 


S 



(5) 


R = 



( 6 ) 


with 


T 


R 



(7) 
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the reduced temperature, and 


M t 


m. m 

i n 

nrT +~m 

i n 


( 8 ) 


the reduced mass; T i} T n , m., m n are the ion temperature, neutral gas tempera- 
ture, ion mass, and neutral particle mass, respectively. Q(g) is given by Stubbe 
(1966) as 


Q(g) 

= 0 

for OlglVj 

Min , 



l V A 


► ^ A -T v o 

(9a) 

Q(g) 


for Vj £ g 00 

J 




Q(g) = 

Q(g) = Q ( 


for 0<g<v„ 


for v 2 — g — 00 


■•(-a 

where v Q , v x , v 2 , Q 0 , and k 0 are defined by 


for A< T v # * 


(9b) 
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CT 2 



!L 

in 


ilk 

f ^ in 


V 

LL . 0 

r in 


( 10 ) 


(ii) 


V — V 
2 0 



Q 0 = 


7T cr* 


V 


k o = 2 % v o - 2ne V~ n 


( 12 ) 

(13) 

(14) 
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e is the ion charge, 77 the neutral particle polarizability, and o- the sum of the 
gas-kinetic radii of the ions and the neutral particles. For A = 0 and a = 0, 
the cross section of (9a) and (9b) becomes identical to the orbiting cross section 
in Langevin's ion mobility theory for an attractive inverse fourth power potential 
(Langevin, 1905). The orbiting cross section was taken as the chemical reaction 
cross section by Gioumousis and Stevenson (1958) in their theory. Introducing 
the further abbreviations 


u 


0 


VR v 2 


(15) 


and 


A 


0 




(16) 


lumping together Equations (3), (4), (9a), (9b) and integrating over g, we obtain 
for the reaction rate constant 


k 


r 



(4Rv s 2 ) n_1 n! 
(2n - 1)! 



v=0 



n=0 


(4Rv s 2 ) n n! 
(2n + 1)! 


£ 

u = n 



V\ 



for 


A > A 0 


(17a) 
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for A < A Q 


(17b) 


For A = A 0 , (17a) and (17b) become identical. When the relative flow velocity 
v s is negligibly small, (17a) and (17b) reduce to 





for A>A 0 


(18a) 


k 



for A<A 0 


(18b) 


-So far we have considered only the conversion of kinetic energy into activation 
energy. There is evidence, however, that energy belonging to the internal de- 
grees of freedom can also contribute to the activation of molecules (e.g. Hinshel- 
wood, 1955; Kontrat’ev, 1964). According to recent studies (Thomas and Norton, 
1966; Schmeltekopf et al., 1967), especially the effect of a vibrational energy 
contribution should be important, and we will, therefore, include it in our 
derivation. 
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The number of molecules being in the j-th vibrational state for a gas in 
thermal equilibrium is given by 

N . = N o (l - e ' h|//1tT ") e “ i ( hl,/kT n) (1 , 

with v the vibration frequency and jhv the energy of the j-th state. Thus, if in- 
stead of the kinetic energy alone the sum of kinetic and vibrational energy is 
available for an activation of the molecule, the reaction rate k r is given by 

kr . •)£ (M <* 

j=o 

k r . (A.) is the reaction rate after (17a), (17b) or (18a), (18b) calculated for an 
effective activation energy 


A j - Max (A - j hv , 0 ) (21) 

instead of A. The vibration frequencies are 4.74 10 13 sec -1 for 0 2 and 7.08 
10 13 sec -1 for N 2 (Syrkin and Dyatkina, 1964). 

We have not yet specified the steric factor, in particular not its tempera- 
ture dependence. In order to do so, we have to leave our collision theory con- 
siderations and engage the more elaborate transition-state theory. The steric 
factor P, as predicted by the transition-state theory, depends on the configuration 
of the activated complex which is formed at the moment the reaction occurs 
(e.g. Frost and Pearson, 1963). We may expect that for both reactions (1) and 
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(2) the activated complex is nonlinear since no valence electrons in the s-state 
are involved. For a nonlinear complex, P is given by 


„ T 7/2 /„ -hv/kT 

3 c \Lli__ 

2 ( 22 ) 

with T c the temperature of the activated complex. A reliable determination of 
T c can be made only when T i and T n are equal. In this case we have, of course, 
T c = T i = T n . Thus, we obtain for the steric factor 

p = p -r-l/2 ^ 

1 1 - e' hv/kT (23) 


0 T. 3/2 T 5/2 


i- 


-hv/WT 


A combination of Equations (20) and (23) with either (17a), (17b) or (18a), (18b) 
gives the final result for k r as a function of the temperature T and, when Equa- 
tions (17a), (17b) are employed, of the relative flow velocity v s . 

The unknown parameters A and P 0 can be determined only empirically by 
selecting them so that the agreement between the theoretical and the experimental 
results is optimum. 

It must be emphasized that the theory developed so far holds true only for 
a system in thermal equilibrium. An attempt to generalize it in order to describe 
the effect of excessively excited vibrational states (Thomas and Norton, 1966; 
Schmeltekopf et al., 1967) will be made in section 5. An attempt will not be 
made, however, to extend the theory to the case T. / T n , because this would 
create many complex problems concerning the partition function and the tempera- 
ture of the activated complex. Another limitation of our theory is set by the 
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fact that for energies above about 0.2 eV (see Hasted, 1964, Chapter 14) some 
of our assumptions will begin to fail. These assumptions are: Applicability of 
classical mechanics, constant gas kinetic radius, pure inverse fourth power 
interaction potential, and well defined geometry of the activated complex. 

The unfortunate consequence of these restrictions is that we cannot make 
use of several valuable experimental works (Bohme, Ong, and Hasted, 1967; 
Giese, 1966; Paulson, Dale, and Murad, 1967; Stebbings, Turner, and Rutherford, 
1966), at least not quantitatively. We can learn from these works, however, that 
despite the possible decrease of k x and k 2 with temperature for low energies 
there is a distinct increase at higher energies. Therefore, our theory must be 
capable of producing a minimum for k f at temperatures in the vicinity of 1000°K. 

3. THE REACTION 0 + + 0 2 - 0 2 + + O 

In this and in the following sections we will present numerical evaluations 
of k r as a function of temperatime. For reaction (1), the polarizability 17 and the 
gas kinetic radius a are taken as 1.60 10 “ 24 cm 3 and 2.50 10 ~ 8 cm, respec- 
tively. Thus, the energy A 0 is 4.72 10 " 13 erg or 6.80 kcal/mole. 

Figure 1 shows k x as a function of T for four activation energies between 
0.86A o and 0.98A o together with the experimental results of Nakshbandi and 
Hasted (1967) and Smith and Fouracre (1968). The data point for T = 375 °K in 

Table 1 of Nakshbandi and Hasted (1967) is disregarded since, as stated by the 

* 

authors themselves, it may be subject to an error. At the time Figure 1 was 
prepared, the author had no knowledge of the work of Dunkin et al. (1968) so that 
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their data points are not included in Figure 1. A redrawing, however, was not 
deemed necessary because of the good agreement between all three experimental 
works. 

The theoretical curves are adjusted to a value of kj = 1.9 10“ 11 cm 3 sec -1 
for T = 300 °K through a proper choice of P 0 . The best over-all fit is afforded 
by the 0.94 curve. The extension towards higher temperatures is presented in 
Figure 2 (solid lines). The 0.94 curve has a minimum at about 650°K. For tem- 
peratures above 750°K, it can be approximately described by 


k, (T) = 0. 58 kj (300) (l + (24) 

The computation of the solid curves in Figure 2 was based on Equation (23) for 
the steric factor P which means that all degrees of freedom of the activated 




Figure 1 Figure 2 
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complex are assumed to be excited up to their full thermal equilibrium extent. 
While this is surely the case for laboratory conditions, it must be doubted that 
in the ionosphere the complex experiences a sufficient number of collisions dur- 
ing its relatively short lifetime to guarantee the full excitation of all degrees of 
freedom, in particular the vibrational degrees. A complete quantitative study of 
this question appears extremely difficult, and it shall therefore not be under- 
taken here. We can, however, roughly estimate the importance of this effect by 
assuming that only one of the two vibrational degrees of freedom is. excited. 

This leads to the dashed curves in Figure 2. 

The assumption that underlies the calculation of the dashed curves in Fig- 
ure 2 is merely based upon the fact that one of the two vibrational degrees of 
freedom, stemming from the molecular reactant, is already fully excited before 
the complex is formed, while the other has to be excited during the lifetime of 
the complex. It is not intended to claim that the situation in the ionosphere can 
be quantitatively described in this way. The only aim of these considerations is 
.to exhibit the possibility of a discrepancy between laboratory and ionospheric 
reaction rates and to estimate its order of magnitude. 

So far all computations were carried out for v s = 0. as a function of v s , 
related to for v s = 0, is shown in Figure 3 (left half) for different gas tem- 
peratures T. We see that k t increases with v s , preferentially for low tempera- 
tures, but that the effect is so small that it should not contribute to any ionos- 
pheric phenomenon. The same holds true for all activation energies below A 0 . 
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However, for activation energies above A 0 , the effect becomes increasingly 
stronger. This is illustrated in the right half of Figure 3 for a hypothetical acti- 
vation energy A = 1.20 A 0 . 

4. THE REACTION 0 + + N 2 -* NO + + N 

As mentioned before, the existing laboratory results for k 2 as a function of 
T are incompatible with each other. We must, therefore, refrain from trying to 
determine A of this reaction and from extrapolating the experimental values for 
k 2 to higher temperatures. What we can do, however, is to give a representation 
of k 2 as a function of both T and A (Figure 4). The results shown in Figure 4 were 
obtained for the parameters a- = 2.55 10' 8 cm, r] = 1.76 1CT 24 cm 3 , and 
A 0 = 4.80 10~ 13 erg or 6.91 kcal/mole. 

The data points of Nakshbandi and Hasted (1967) can approximately be re- 
produced by a curve lying between curves 5 and 6 of Figure 4. This means that 
k 2 increases monotonically with T and reaches a value at 2000°K 1.65 times that 
at 300 °K. 
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Figure 4 


The results of Dunkin et al. (1968), on the other hand, can be approximated 
by a curve in the vicinity of curve 1 in Figure 4. Correspondingly, k 2 has a weak 
minimum around 1100°K and retains a nearly constant value of about half that at 
300°K up to 2000°K. Thus, a reaction rate of k 2 % 1.0 10~ 12 cm 3 sec -1 for 
T = 300°K (Dunkin et al., 1968) should result in a value of k 2 between 0.5 and 
0.7 10" 12 cm 3 sec -1 for F-region temperatures. 


5. THE EFFECT OF ENHANCED VIBRATIONAL EXCITATION OF N 2 ON 
THE REACTION 0 + + N 2 - NO + + N 

The theory presented in section 2 is bound to the restrictions Tj = T v (T v 
= vibrational temperature of N 2 ). Therefore, modifications are necessary in 
order to use it for a description of the experimental results of Schmeltekopf 
et al. (1967) who find that for a gas temperature T = T. = T n = 300 °K the reac- 
tion rate k 2 experiences an increase by a factor of 30 when T v is changed from 


300 °K to 6000°K. 
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Equation (20) can easily be generalized to 


k 


r 


i=o 


-j(hv/kT v ) 


k rj ( A r T ) 


(25) 


with A i and k r . given by Equations (21) and (17a), (17b) or (18a), (18b), respec- 
tively. The difficulty lies in the steric factor, i.e. in the temperature T c of the 
activated complex. After (22) , P is given by 

i - .- h ” /kT ') 


T m I 

P = p 0 1 


(l ’ 


(26) 


A straightforward purely theoretical determination of T c appears extremely 
intricate. However, T c can be readily estimated in the following semi-empirical 
way: For T = T v , the average internal energy of the activated complex is 


U c = \ kT + 2hv (l - e" hlVkT ) ^ j e _j ( hv/kT > 

j = i 


(27) 


If T v is greater than T, the activated complex is supplied with a fraction a of the 
excessive vibration energy of N 2 , 


All 


= h, (l - e"”""’ ) J2 j - hv(l - e*/«) £ 


J e 


-Xhv/kT) 


(28) 


Thus, T c can be obtained from the determining equation 


u c * 


j=l 


(29) 


a has to be determined empirically by selecting it so that the theoretical result 
gives the best fit to the experimental result. 

Figure 5 shows k 2 as a function of T v . The thicklined curve represents the 
smoothed data of Schmeltekopf et al. (1967). The thinlined curves are the theo- 
retical results for four gas temperatures between 300°K and 2000°K, computed 
for a = 0.40 and A = A 0 . It is clearly recognizable that the increase of k 2 with 
T v is much more pronounced for low temperatures than it is for higher tempera- 
tures. This means that the effect of an enhanced vibrational excitation of N 2 on 
ionospheric F-region phenomena should be much stronger for nighttime and solar 
minimum conditions than for daytime and solar maximum conditions, respectively. 



Figure 5 
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The choice of the activation energy A = A 0 isto some extent arbitrary. 
Computations for different activation energies were carried out. They showed 
that a slight dependency on A exists but that in no case the qualitative character- 
istics, as described aboye, were changed. Thus, we can and will put emphasis 
only on these qualitative features. 

6. CONCLUSIONS 

In the present paper an attempt is made to theoretically describe the reac- 
tion rate k r of an ion -molecule interchange reaction as a function .of the tempera- 
ture T, the relative flow velocity v s of the gases, and the parameters A (activation 
energy) and P Q (constant part of the steric factor P). The resulting formulae are 
applied to the reactions 0 + + 0 2 - 0 2 + + O and 0 + + N 2 - NO + + N. 

For the first reaction, the rate constant is found to have a minimum at 
about 650°K and a nearly linear increase for the range of ionospheric tempera- 
tures. The possibility of a discrepancy between laboratory and ionospheric 
reaction rates is discussed. It is shown that kj depends on v s only weakly but 
that the dependency can become important for higher activation energies. 

For the second reaction, the experimental situation is contradictory. There- 
fore, a determination of A becomes impossible, and k 2 is presented as a function 
of both T and A. For further experiments it may be suggested to emphasize low 
temperature measurements since they contain more information about this reac- 
tion than high temperature data (see Figure 4) . 
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By introducing an additional parameter, the theory is extended in order to 
make it applicable for a study of the effect of enhanced vibrational excitation of 
N 2 on the reaction rate k 2 . It is found that the effect strongly decreases with 
increasing gas temperature T, Therefore, it should be more important for the 
ionosphere under nighttime and solar minimum conditions than under daytime 
and solar maximum conditions. 
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FIGURE CAPTIONS 


Figure 1. k % as a function of T for the activation energies A/A 0 = 0.86, 0.90, 

0.94, 0.98. Triangles: Experimental results of Nakshbandi and 
Hasted (1967). Circles: Experimental results of Smith and Fouracre 
(1968) 

Figure 2. Ibid, for Figure 1. Solid curves: Both vibrational degrees of freedom 
of the activated complex excited. Dashed curves: Only one vibrational 
degree of freedom of the activated complex excited 

Figure 3. k x as a function of v s . Left half for A = 0.94A 0 , right half for 
A = 1.20A 0 

Figure 4. k 2 as a function of T for nine different activation energies between 
0.80A 0 and 1.20A o 

Figure 5. k 2 as a function of the vibrational temperature of N 2 , T v . Thick- 
lined curve: Experimental results after Schmeltekopf et al. (1967) 
for a gas temperature T = 300°K. Thinlined curves: Theoretical 


results for the gas temperatures 300, 500, 1000, 2000°K 



